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<5> Abstract 

o 

We consider density functionals for exchange and correlation energies in two-dimensional systems. The functionals are constructed 
by making use of exact constraints for the angular averages of the corresponding exchange and correlation holes, respectively, and 
^ assuming proportionality between their characteristic sizes. The electron current and spin are explicitly taken into account, so that 
O the resulting functionals are suitable to deal with systems exhibiting orbital currents and/or spin polarization. Our numerical results 
show that in finite systems the proposed functionals outperform the standard two-dimensional local spin-density approximation, 
still performing well also in the important limit of the homogeneous two-dimensional electron gas. 
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X. Introduction 

' Practical success of density-functional theory fl*] (DFT) de- 
pends on finding good approximations for exchange and cor- 
relation energy functionals. So far, most density functionals 
have been developed in three dimensions (3D) with a view to 
studying the properties of atoms, molecules, and solids. Apart 
from very recent progress 12. B 0. S S S Si. such efforts 
for two-dimensional (2D) systems have been relatively scarce 
beyond the commonly-used 2D local spin-density approxima- 
tion (LSDA). However, the rapidly increasing theoretical, ex- 
perimental, and technological interest in 2D structures such 
as semiconductor layers and surfaces, quantum Hall systems, 
graphene, and various types of quantum dots, calls for further 
developments. 

Here we consider expressions of density functionals for ex- 
change and correlation energies of 2D systems. These approx- 
imations are derived |l2|, |3t] by modeling the angular average 
of the exchange and correlation holes and using first-principle 
arguments along the lines of Refs. ||^[l3l- The resulting spin- 
and current-dependent density functionals are then combined to 
obtain a new approximation for the exchange plus correlation 
energy beyond the LSDA. Numerical tests, performed on sys- 
tems for which accurate reference data is available, show that 
our functionals are superior to the LSDA. Moreover, we find a 
good agreement with the exact result for the homogeneous 2D 
electron gas (2DEG). 

2. Basic formalism 



Within the Kohn-Sham (KS) method of spin-DFT, UM the 
ground state energy and spin densities p-f(r) and px(r) of a sys- 
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tem of N — N-\ + Ni interacting electrons are determined. The 
formalism is general in terms of dimensionality and the form 
of interaction. However, here we focus on Coulomb-interacting 
particles in 2D. The total energy, which is minimized to obtain 
the ground-state energy, is written as a functional of the spin 
densities (in Hartree atomic units) 



^ I drv^(r)pAr), 

,--1 I *J 



(1) 



o-=T,l 



where Ts\p'\,p[\ is the KS kinetic energy functional, Vo-(r) is 
the external (local) scalar potential acting upon the interacting 
system, £'h[p] is the classical electrostatic or Hartree energy of 
the total charge density p(r) - p-r(r) -i- pj,(r), and Exc\p-\,Pi\ is 
the exchange-correlation energy functional. The latter can be 
further decomposed into the exchange and correlation parts as 



Exc\P],Pi\ = ^^[PT.Pi] + £'c[PT.Pl]- 



(2) 



This is the central quantity which is formally well defined and 
exact, but needs to be approximated in practical applications of 
DFT. 



3. Modeling the exchange 

The exchange-energy functional can be expressed as 



Ex\P],Pi\ 



where 



r 

Jo 



t/^(r) = 27r dsh°;(r,s) 



(3) 



(4) 
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is the exchange-hole potential, and 

-In 



(5) 



is the 2D cylindrical average of the exchange-hole function 



/2:^(ri,r2) = - 



(6) 



where the sum in the numerator is the one-body spin-density 
matrix of the (ground state) Slater determinant constructed from 
the KS orbitals, ipk.a- As the basis of our exchange model 
we consider 



h'^^ia, b; s) = - - exp [-a {b + s^')\ Mia V^s) , 



(7) 



where Iq{x) is the zeroth order modified Bessel function of the 
first kind (note that /o(0) = 1). This function provides a non- 
positive exchange hole and satisfies the normalization 
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In Eq. (|7]i, a and b are position-dependent functions, which 
are chosen to locally reproduce the short-range behavior of the 
exact exchange hole. Comparison to the Taylor expansion leads 
to relations a - npa-e^' and b - y/a, where y :- ab is the 
solution of 



(y- l)exp(y) = 



where 



1 C^Paf 



2 Pa 



+ 2 



(9) 



(10) 



is the so-called curvature of the exchange hole, Tq- is (twice) 
the spin-dependent kinetic-energy density, and o- is the spin- 
dependent paramagnetic current density. Both t,j- and j;,,,,- de- 
pend explicitly on the KS orbitals. Thus the expression in 
Eq. ( [Tol l has an implicit dependence on the spin-densities p,j- 

Note that the model introduced in Eq. O is exact in the case 
of the ground-state single-electron wave function of a 2D har- 
monic oscillator Furthermore, application of the above model 
to the 2DEG yields the exchange energy per particle as 



ejr,,^] 



4V2 



[(1+^)3/2 + (l_^)3/2j^ (11) 



where r, = 1/ s/jrp and ^ - (pi - pi)/p is the spin polarization. 
This energy is about 4.4% smaller than the exact one lIlSIl . Inter- 
estingly, Eq. (fTTT i may be seen as a correction of the LSDA ex- 
change functional suitable for dealing with finite few-electron 
systems. 

4. Modeling the correlation 

The correlation-energy functional can be expressed as fol- 
lows 

E,[p^,Pi] = El^[p^,pi] + £iJ-[pT,Pi] + 2El^[p^,pi], (12) 



T I IT 

where we have used -Ec [Pt'PJ.] - [Pt'PJ 

E^^'lpj,Pl]^n jdrpjr) dsh-'/{r,s) (13) 

with /!^°" (r, s) being the cylindrical average of the correlation- 
hole function, /i^"" (ri, r2), computed in the same way as 
Jf^iv, s). We remind that the correlation-hole functions can be 
obtained by the coupling constant integration 



(14) 



where /j^^ (^^^2) = /z^ (ri,r2) - /ij(ri,r2)5<,^s {crcr'} = {TT 
, ii, Ti, it)^ and the parameter A e [0, 1] is the electronic cou- 
pling strength 1 1]. Now, we model [3] the cylindrical averages 
of the /l-dependent correlation-hole functions as 



rX(r,i)= -^D^(r) 



and 



hf,{r,s)^2Apdr) 



1 + iAzaAr) 



s^F{yaAr)s), (15) 



1 + lAzaAr) 



(16) 



for the same spin and opposite spin case, respectively (note 
that, (T indicates an opposite value w.r.t. cr). In Eq. (fTST l. 
Da- :- 1/8 (V^po- - 4Cj). This quantity vanishes for all the 
single-particle (A^ = 1) systems, and thus also the correspond- 
ing (same-spin) correlation hole and correlation energy vanish 
as in the exact case. For 5 — » 0, Eqs. (fTsT i and ( fTSI l recover the 
exact short-range behavior of the /l-dependent correlation-hole 
functions. F is introduced to ensure the decay of the correlation 
holes in the limit i — > 00. In the case of finite 2D systems, it 
seems appropriate to set F(x) :- exp(-x^). Furthermore, the 
exact condition 

(17) 



is satisfied by setting 7^0- = 3 y/n/iAzaa) andjaa = y/n/ilzaa-)- 
Functions Zaa and Zaa are defined as the characteristic sizes of 
the correlation holes, for which /j^^(r, Zo-a) = 0, h'^'^(r, Zaa-) = 
0. As for the 3D case, 10 we set Zaa IcaaW^f^ and 
z^g. Cct-o- [|f/.?"r' -i-|f/fr'] where f/^ is the exchange-hole 
potential calculated in the previous section, and Caa and Ca-a 
are constants to be determined (see below). In this way, the 
characteristic size of the correlation hole is proportional to the 
characteristic size of the corresponding exchange hole. This is 
suggested by a picture in which electrons are correlated as long 
as their interaction is not screened by their respective exchange 
holes. A more elaborate model may require to consider Caa' 
as a density functional, which is, however, beyond the scope of 
this work. 

The approximations given above are implicit density func- 
tionals due to the explicit dependence on the KS orbitals ap- 
pearing in the expressions for t and jp o-. Therefore, these func- 
tionals can be classified as current-dependent meta-generalized- 
gradient approximations, and we expect them to be ideally 
suited to deal with a large class of spin-polarized and/or current- 
carrying ground-states of two-dimensional systems. 



5. Applications 

We test the present 2D functionals for a set of A'-electron 
parabolic quantum dots described by a Hamiltonian, 
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(18) 



where oj - 1 /4 is the confinement strength. The exchange en- 
ergies obtained using the approximation given in Sec. |3] are 
compared with the exact-exchange (EXX) resuhs within the 
Krieger-Li-Iafrate approximation |12]. Then, the correlation 
energies obtained according to Sec. |4] are compared to the ref- 
erence results El^* = £[^1 - Ef^^, where Ef^f is the numerically 
accurate total energy of a configuration-interaction calculation 
reported in Ref. [14]. TheKS orbitals from self-consistent EXX 
calculations - performed using the octopus DFT code lll3ll - 
are used as the input for the functionals. 

In Table[T]we show the results for fully spin-polarized quan- 
tum dots with different = 2 ... 6, spin S - N/2, and total 
angular momentum (z component) L. Note that different val- 
ues of L correspond to different paramagnetic current densities. 
Overall, we find excellent performance in both the exchange 
and correlation functionals. A significant improvement over the 
LSDA as seen in the mean errors A on the last row of Table[T] It 
should be noted that the error of the LSDA correlation is huge. 
This error is typically partially compensated by an error of an 
opposite sign in the LSDA exchange. However, our functionals 
are able to reproduce both the exchange and correlation ener- 
gies with a good accuracy. 

In the correlation functional we have used a fixed value 
Co-cr - 1 -25 as the prefactor in the characteristic size of the cor- 
relation hole (see the end of Sec.|4]i. We point out that although 
this value is fundamentally a post-hoc choice, the variations in 
optimal Cg-cr (and Cg-g-) for different systems are typically rather 
small 1 3]. This becomes obvious also within the 2DEG below, 
where results for various densities, corresponding to amounts 
of correlation, are shown with the same fixed Cg-a- 

Finally we turn our attention from finite systems to the 
2DEG. Figure \T\ shows the exchange (left) and correlation en- 
ergies (right) per particle in the 2DEG as a function of the den- 
sity parameter r, = 1 / against exact results available from 
analytic |15] and quantum Monte Carlo calculations |16], re- 
spectively. As pointed out in Sec. [3] in this limit our exchange 
functional deviates consistently 4.4% from the exact (LSDA) 
result. Also for the correlation the agreement is very good up to 
strong correlations - it should be noted that in quantum dots, for 
example, typical densities correspond roughly to 1 < ^ 5. In 
the correlation functional we have consistently used fixed val- 
ues Ca-cr - 1.25, i.e., the same value as with the quantum dots 
above, and Co-o- - 0.75. 



6. Conclusions 

We have presented two-dimensional density functionals for 
exchange and correlation energies. The functionals have been 
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Figure 1: (color online) Exchange (left) and correlation energies (right) per 
particle as a function of the density parameter r, = 1 / V™ in the fully spin- 
polarized and spin-unpolarized homogeneous two-dimensional electron gas. 
The dashed lines show the results calculated with the present functionals. The 
solid lines con'espond to the exact results obtained from analytic calculations 
for the exchange |15] and analytic parametrizations of the quantum Monte 
Carlo data for the correlation 1 16]. Note the logarithmic scale in the exchange 
energies. 



derived by modeling the angular averages of the correspond- 
ing exchange and correlation holes along first-principle argu- 
ments and by enforcing known exact constraints. Tests of our 
functionals for few-electron open-shell quantum dots and the 
homogeneous two-dimensional electron gas show very good 
performance. Hence, we may propose an accurate alterna- 
tive to the commonly used local spin-density approximation for 
electronic-structure calculations in two dimensions. 
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Table 1: Exchange and correlation energies for fully spin-polarized quantum dots. The last row shows the mean percentage error A of out functional and the LSDA 
for both exchange and correlation, respectively. 
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